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ABSTRACT 

We extend recent numerical results ([Dobbs et al.ll2QQ6l ) on molecular cloud formation 
in spiral galaxies by including a multi-phase medium. The addition of a hot phase of 
gas enhances the structure in the cold gas, and significantly increases the fraction of 
molecular hydrogen that is formed when the cold gas passes through a spiral shock. 
The difference in structure is reflected in the mass power spectrum of the molecular 
clouds, which is steeper for the multi-phase calculations. The increase in molecular 
gas occurs as the addition of a hot phase leads to higher densities in the cold gas. In 
particular, cold gas is confined in clumps between the spiral arms and retains a higher 
molecular fraction. Unlike the single phase results, molecular clouds are present in the 
inter-arm regions for the multi-phase medium. However the density of the inter-arm 
molecular hydrogen is generally below that which can be reliably determined from CO 
measurements. We therefore predict that for a multi-phase medium, there will be low 
density clouds containing cold atomic and molecular hydrogen, which are potentially 
entering the spiral arms. 

Key words: galaxies: spiral - hydrodynamics - ISM: clouds - ISM: molecules - stars: 
formation 



1 INTRODUCTION 



Numerical simulations have become a powerful tool to ex- 
. plore the physic s of star formation (lOstriker et aL I2OOII: 
Bate et all 12003': 'Mac L ow KlessenI I2OO4": B onnell et all 



20071 ). Recent advances in computational power now make 
it possible to study the formation of giant molecular clouds 
(GMC) and the onset of star formation (iKim et al ] I2OO3I: 
Bonnell et all I2OO6I: IPobbs et all [2006| : iGlover Mac Lowl 



2006al lbl: Vazquez- Semade ni et al.] r2006V There are two gen- 
eral formation mechanisms that have been advanced, the 
first is that a cloud becomes self- gravitating, inducing high- 
densities and the conversion of HI into molecular gas. The 
second possibility is that an external triggering such as a 
spiral or other turbule nt shock either gath ers together pre- 
existing molecular gas (|Pringle et al.|[200l[ ) or induces high 
densities in the shock and thus converting HI into H2. 

Previous results show that spiral shocks induce the for- 
mation of dense molecular clouds (Dobbs et al. 2006) and si- 
multaneously generate a velocity dispersion (Bonnel l et al.l 
1 2006 ). A prerequisite for the formation of these molecular 
clouds in isothermal simulations is that the gas is cold, of 
order 100 K or less (IPobbs et al.ll2006l ). For cold gas, the 
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densities in the shocks are sufficiently high for H2 forma- 
tion, and the dynamics of the spiral shocks lead to clumpy 
structures in the gas. The cold atomic component of the 
ISM contains the most mass of atomic hydrogen in the ISM, 
but the fi lling factor of this phase is very small 0.013, 
[C^l2005l ). The warm (lO^K) atomic phase, which is typi- 
cally used in simulations of spiral galax ies (Ki m Ostrikerl 
l2002l : IWada k Kodal[20o3 : iGittins Clarkell20o3 ) . contains 
less mass but f ills just over a third of the volume of the 
ISM (|Coxll2005l V The simulation described in this paper is 
the first to examine the dynamics of a spiral potential with 
both cold and warm phases. The resulting gas distribution 
is more representative of the ISM, with cold gas embedded 
in smaller clumps within a warm diffuse phase. 

In the present study, we investigate the effects of the 
hot medium on a pre-existing cold HI as it passes through 
the spiral shock. The motivation is to include the pres- 
sure confinement of the hot gas in the dynamics of a spiral 
shock. We do not include any transition between the hot 
and cold c omponents but rath er assume the picture pre- 
sented by (jFringle et al.1 l200lh that molecular clouds are 
formed from gas that is cold before it enters the shock. 
We leave a full treatm ent of the thermal e volutio n of the 
gas to a later study. IGlover Mac LoJ (|2006al lbl) ana- 
lyze molecular cloud formation over much smaller size- 
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scales (20-40 pc), but include a more consistent treatment 
of the ISM and H2 formation. In their models, molecular 
clouds form slowl y in gravitationally collapsing uniform gas 
([Glover Mac L ow 2006 j) or rapidly from t he compression 
of turbulent gas (Glover Mac Lo w 2006b). The gas tem- 
perature in these calculations is initially 10"^ K, but quickly 
drops to 100 K or less once the density exceeds 1 cm~^. 
Other models of the ISM have shown the formation of cold 
(< 100 K) HI, as a precursor t o molecular clouds, on local 
(IVa zguez-Semad eni et all Il995l : |Audit & Hennebelle 2005; 
[Hdtsch et al. 20Q5h and galactic ( Wada Norman 1999.) 
scales. 

We present numerical results to investigate the dynam- 
ics of cold gas subject to a spiral shock when in the pres- 
ence of a confining warm phase. We compare this 2 phase 
simulation of a galact ic disc with t he pr evious single phase 
simulations shown in llDobbs et al.l (|2006h . In particular, we 
describe the structure of molecular clouds and the content of 
molecular hydrogen across the galactic disc. We also provide 
the first attempts to determine individual molecular cloud 
properties from global simulations. 



2 CALCULATIONS 

We use the 3D smoothed partic le hydrody namics (SPH) 
code based on the version by Benz (|Benzlll990l l. The smooth- 
ing length is allowed to vary with space and time, with the 
constraint that the typical number of neighbours for each 
particle is kept near Nmigh ~ 50. Artificial viscosity is in- 
cluded with the standard parameters — 1 and [3 — 2 
(|Monaghan k Lattanz'iolll985l : lMonaghanlll992l ). 



2.1 Flow through galactic potential 

The galactic potential includes spiral components for the 
disc, dark matter halo and spi ral density pattern . The 
potential is described in full in iDobbs et al.[ (|2006l ). The 
spiral component h as 4 spiral arms and is taken from 
I Cox &: GomezI (|2002h . The amplitude is 1 atom cm~^, the 
pattern speed is 2 x 10~^ rad yr~^ and pitch angle a — 15^^. 
The pattern speed leads to a co-rotation radius of llkpc. 
The disc is in equilibrium, as the rotational velocities of the 
gas balance the centrifugal force of the potential. 

Self-gravity magnetic fields, heating, cooling or feedback 
from star formation are not included in these calculations. 
We further neglect pressure due to cosmic rays, which con- 
tributes tot the pressure of the diffuse ISM. The paper in- 
stead focuses on how hydrodynamic forces and the galactic 
potential influence the flow. 



2.2 Initial conditions and calculation of molecular 
gas density 

Gas particles are initially distributed as described in 
iDobbs et all (|2006h . occupying a region of radius 5 kpc< 
r <10 kpc. The disc has a scale height z ^ 100 pc, and a 
velocity dispersion of 2.5% of the orbital speed for the z com- 
ponent of the velocities maintains vertical equilibrium. We 
determine the x and y components of the positions and ve- 
locities from a 2D test particle run. Particles in the test run 
are initially distributed uniformly with circular velocities. 



They evolve for a couple of orbits subject to the galactic po- 
tential, to give a spiral density pattern with particles settled 
into their perturbed orbits. A velocity dispersion of 2.5% of 
the orbital speed is also added to the x and y velocities. 

The gas is distributed in 2 phases, a cold component 
of T = 100 K, and a warm component of T = 10^ K. 
Each phase constitutes half of the gaseous mass of the disc, 
5 X 10^ M0, so the total mass of the disc is 10^ M©. This 
corresponds to a nominal surface density of ~ 4 M© pc^ 
(9 X 10~^ g cm~^). The total HI surf ace density at the solar 
radius is - 5 M© pc^ (Wolflre et a ll [2003 ). The SPH parti- 
cles are initially randomly assigned as hot or cold gas, but 
the distribution settles into a diffuse warm and dense cold 
phase as the simulation progresses. The gas is isothermal 
for both phases. The calculation of the smoothing lengths 
does not distinguish between cold and hot particles. Using 
two phases can lead to numeric al issues with the smooth- 
ing lengths (|Thacker et al ]l2000'), in particular overcooling, 
although both phases are isothermal in our simulations. Fur- 
thermore only a small fraction of cold particles have predom- 
inantly hot neighbours e.g. only 2 % of cold particles have 
< 20% cold neighbours. 

Ideally, a full treatment of the ISM, which includes heat- 
ing and cooling between the two phases should be applied. 
However, the primary aim of this paper is to investigate 
how a hot phase affects the dynamics of cold gas passing 
through a spiral shock. There have been recent theoretical 
arguments that molecu lar clouds ma y form from cold HI 
or molecular gas (Prin gle et al.|[200lh , motivated by obser- 
vational evidence for shorter cloud lifetimes. Observations 
also i mply that m olecular clouds form from molecular gas in 
M51 (jVogel et al.l ll988). and thre may be reservoirs of low 
density H2 present in galaxies which is difficult to detect 
(iLeaueux et al. 1993: Allen et al. 2004). Cold HI is ubiqui- 
tous in the Outer Galaxy ( Gibson et al.ll2006ll and ap pears 



tous m tne Uuter ualaxy (lljibson et aLllzUUbl) and ap i: 
to be flowing into the spiral arms ([Gibson et al ] l2005b . 

We calculate the density of molecular hydro gen by the 
same process as described in iDobbs et all (|2006h . We post- 
process our results to determine the fraction of molecular 
gas, using the equ a tion f or the rate of change of H2 density 
from lBergin et al.1 (|2004l ) 



dn[H2 



Rgr{T)npn{H) - [Ccr + C,diss{N{H2),Av)]n{H2 



(1) 

The total number density is rip, where Up — n(H) +2n(if2), 
N is the column density (of atomic or molecular hydro- 
gen) and T the temperature. The strength of the UV fleld 
is Cdiss{0) = 4.17 X 10~^^ s~\ comparable with the UV 
flux from a BO star, and the cosmic ray ionization rate is 
= 6 X 10"^^ s"^ There is some uncertainty in the cos- 
mic ray ionization rate (e.g. iLe Petit et al.l 120041 1 but this 
term is much less than the the other components of Eqn. 1 
and unUkely to affect our results. Rgr is the formation rate 
on grains, assuming an efficiency of S = 0.3. We do not ex- 
plicitly change S for the warm gas in our calculations, but 
the densi ty of the wa rm gas is too low for H2 formation 
to occur (Pobbs et al.|[2006 ) and H2 is dissociated at high 
temperatures. Ideally, a more realistic treatment of the ISM 
including heating and cooling should be included for a more 
ac curate calcul a tion of H2 formation. However, as shown 
m toobbs et al. this method primarily extracts the 

densest gas. 
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Our main simulation was performed with 6 million par- 
ticles, 3 million in each of the hot and cold phases. This 
compares to o ur previous sing:le-p hase simulation of 4 million 
cold particles ([Pobbs et al.|[2006l ). In addition to these high 
resolution runs, we ran lower resolution runs (of 1 million 
and 250,000 particles) to investigate the effects of the nu- 
merical resolution. The lower resolution runs generally yield 
lower maximum densities throughout and lower H2 fractions. 



3 RESULTS 

We first discuss the overall structure of the disc and the 
overall distribution of molecular hydrogen. In Section 3.5, 
we use a clump-finding algorithm to determine the prop- 
erties of individual molecular clouds. We compare results 
from the multiphase calculati on and t he sin gle phase (50 K) 
simulation described in Dobb s^t al.l (|2006l ). In Sections 3.3 
and 3.6.1 we compare the inter- arm and arm molecular gas 
densities for the multi-phase simulation with observational 
results, and comment on the ratio of inter-arm to arm molec- 
ular clouds. 

3.1 Structure of the disc 

The column density section of the disc of the multi-phase 
simulation is shown in Fig. 1 at different times. From Fig. 1, 
we see the growth of small scale structures into self consis- 
tent, large scale features over time. The figure includes both 
the 100 and 10^ K gas. The 100 K gas is most clearly visible 
in small clumps whilst the 10^ K gas provides a diffuse back- 
ground and the smooth component of the spiral arms. The 
2 components are shown explicitly in Fig. 2, which plots the 
particles representing the hot and cold gas separately for a 
section of the disc after 100 Myr. Both figures show that 
the cold gas is confined by the pressure from the hot com- 
ponent into small dense clumps. The hot gas on the other 
hand is much more uniform and smooth. As previously de- 
scribed (jPobbs et al.ll2006h . the structure of the gas in the 
spiral arms is dependent on the gas dynamics of the spiral 
shock. For cold gas, the spiral shock increases clumpy struc- 
ture in the gas. For hot gas the shock is much weaker, so 
has much less effect on the dynamics, an d the gas pressure 
will s mooth out any structure in the gas (jPobbs Bonnell 
I2OO6I ). Two further differences are apparent for the spiral 
shock for each component. As can be seen in Fig. 2, the 
shock for the cold gas is offset from the shock corresponding 
to the hot gas. However the shock from the hot component 
is evidently affecting the cold gas clumps, as can be seen 
from the elongation of the cold gas clumps as they enter the 
spiral arm (Fig. 2). The shock for the hot gas is also much 
wider than for the cold gas, since the degree of compression 
of the hot gas is much less, whilst the cold gas forms a very 
narrow spiral arm. 

Fig. 3 shows the distribution of mass with density for 
the multi and single phase simulations after 100 Myr. The 
hot component contributes a greater mass at low densities 
for the multi-phase simulation. The greater mass at high 
densities, compared to the single-phase disc is due to the 
confinement of the cold gas to higher densities by the hot 
component. 
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Figure 1. Column density plots (g cm~^) showing a 5 kpc by 5 
kpc section of the disc (with Cartesian coordinates -9 kpc < x < 
-4 kpc and -5 kpc < y < kpc. The xy-coordinate grid is centred 
on the midpoint of the disc and remains fixed with time. Gas is 
flowing clockwise across the disc, i.e. from bottom to top. The 
time corresponding to each plot is i) 60 Myr, ii) 140 Myr and 
iii) 220 Myr. The plot shows the total column density, including 
both the 100 and 10^ K components. 



Whilst the distribution of hot gas changes little with 
time, the structure of the cold gas varies. Until approxi- 
mately 60 Myr (Fig. li), the cold gas is fairly regularly dis- 
tributed in small clumps across the section of disc, with an 
increased concentration in the spiral arms. Unlike the single- 
phase simulations, where the gas is clumpy due to the parti- 
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Figure 2. The particle positions are shown for a section of the 
disc, assuming a Cartesian grid centred on the midpoint of the 
disc. The 2 phases are plotted, 100 K (blue) and 10"^ K (red) after 
100 Myr. 



cle nature of the code, the warm phase here confines the cold 
gas into dense clumps. For the single-phase calculations, the 
dumpiness in the spiral arms is amplified with time by the 
dynamics of the shock (jPobbs et al ]l20p6l). In these multi- 
phase calculations however, the spiral shock acts rather to 
organise these clumps into more coherent structures. These 
are perhaps more easily visible as the larger features shear- 
ing away from the spiral arms in Fig. lii and iii. By 220 Myr 
(Fig. liii) the 100 K spiral arm clumps have been sheared 
into thin filaments between the arms, whilst new structures 
are breaking away from the spiral arm. 



3.2 Molecular gas density 

We show column density plots for the simulation after 100 
Myr in Fig. 4. The column density of molecular hydrogen 
is highlighted in red. Fig. 4 shows more clearly the differ- 
ence in structure of the cold gas in the single and multi- 
phase simulations, in particular by comparing Fig. 4 iv) and 
the equivalent frame of the single-phase simulation (Fig. 9 
iv, iDobbs et al.l l2006l l . The multi-phase gas is much more 
structured, with many smaller scale clumps visible. There 
are also small clumps of cold gas strewn across the inter-arm 
regions. Gas in the single-phase int er- arm regions is sheared 
into spurs and feathering (Dobbs Bonnellll2006l ) . but indi- 
vidual clumps and structures are much more defined in the 
multi-phase medium. The cold gas assembles into longer, 
more coherent, very thin spurs in the inter- arm regions at 
later times in the multi-phase calculation (Fig. liii). 

The molecular hydrogen corresponds to 12% of the total 
mass, i.e. 1.2 xlO^ M© and 24% of the cold gas. Most of the 
clumps of cold gas are over-plotted with the molecular hy- 
drogen column density. T he most striking diff erence between 
these results and those of iDobbs et aD (|2006h (Fig. 9) is the 
greater amount of molecular gas. The fraction of molecular 
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Figure 3. The distribution of mass as a function of gas density is 
shown after 100 Myr for the single (crosses) and multi-phase (100 
K, squares and 10^ K gas, triangles) simulations. The nominal 
densit y for molecular cloud s is n(H2)=100 cm""^ or 3 x 10~^^ g 
cm-3 dWilliams et al.ll2000h . 



gas is approximately twice that of the single-phase simula- 
tion after 100 Myr. Unlike the single-phase calculations, the 
cold gas is confined by the pressure from the hot gas and 
remains more dense in the inter-arm regions. Consequently 
there is less photodissociation of H2, and higher densities 
of molecular gas are produced between the spiral arms. As 
seen from Fig. 4 iii) and iv), the molecular gas does not be- 
come fully dissociated between the spiral arms, unlike the 
single-phase medium. The mass of molecular gas continues 
to rise throughout the multi-phase simulation, unlike the 
single-phase results, which peaked after about 100 Myr. As 
the gas is not fully dissociated between the arms, the total 
molecular gas density continues to increase after multiple 
spiral arm passages. By the end of the simulation the frac- 
tion of molecular gas appears to be converging to approxi- 
mately 25 % of the total mass or equivalently 50 % of the 
100 K gas (2.5 x 10^ M©). These figures are however depen- 
dent on the scale height used for the column densities of HI 
and H2, which was assumed to be 100 pc (the scale height 
of the disc). Taking a scale height of 25 pc (comparable to 
the nearest B star from the Sun) reduces the total mass of 
H2 by 30%. 

Fig. 5 shows the density of the 100 K gas (total HI 
and predicted H2), the predicted H2 density and the 10^ K 
gas density as a function of azimuth. These are the average 
densities, calculated over a ring centred at 7.5 kpc, of width 
200 pc and divided into 100 segments. This figure again 
shows the much higher density of H2 in the inte r- arm regions 
compa red to the single-phase medium (Fig. lO jPobbs et al.l 
(|2006 ^)) with the density of H2 now 1/lOth of the total cold 
HI + H2 gas in the inter-arm regions. The mass of molecular 
gas is now approximately 1/3 of the total mass of cold gas 
in the spiral arms, and 1/5 in the inter-arm regions. Fig. 5 
further shows the difference in location of the shocks, as the 
density of the hot gas peaks at a lower azimuth than the 
cold gas. 

Fig. 6 shows the typical timescale gas spends with dif- 
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Figure 4. Column density plots (g cm~^) showing density of molecular hydrogen (red) against overall density (black and white). Length- 
scales of plots are i) 20 kpc x 20 kpc, ii) 10 kpc x 10 kpc, iii) 5 kpc x 5 kpc (-8.5 kpc < x < -3.5 kpc, -8.5 kpc < y < -3.5 kpc), iv) 3 
kpc X 3 kpc (-6.5 kpc < x < -3.5 kpc, -6.5 kpc < y < -3.5 kpc). The xy-coordinates assume a Cartesian grid centred on the midpoint of 
the disc which remains fixed with time. Gas is flowing clockwise across the disc. 



ferent fractions of H2. The average time for each fraction 
of H2 is determined from the duration particles exceed that 
fraction of H2 during each simulation. The typical time gas 
contains over 50% molecular hydrogen is around 40 Myr for 
the single phase simulation and 60 Myr for the multi-phase 
simulation. These times decrease to approximately 5 and 10 
Myr for gas which contains over 90% molecular hydrogen. 
When a scale height of 25 pc is used, these timescales de- 
crease by around 40%. 

3.3 Inter-arm atomic and molecular gas 

The distribution of mass with density is also shown specif- 
ically for the cold component in the inter-arm regions in 
Fig. 7. Clearly the inter-arm gas in the multi-phase simu- 
lation reaches densities over a magnitude higher than for 
the single phase simulation. From local simulations which 
include cooling, gas with densities > 10 cm~^ is found to 



be col d, i.e. 100-200 K or less (Fig. 1, iGlover Mac Lowl 
2006a). We therefore expect a significant proportion of the 
cold gas entering the spiral arms to remain cold in a multi- 
phase medium. By contrast, for the single phase calculation, 
all the gas in the inter-arm regions is likely to increase in 
temperature above 100 K. 

For the multi-phase simulation, the density of the lOOK 
gas (inclusive of H2) is on average 1 — 2 x 10~^^ g cm~^ and 
7 X 10~^^ g cm~^ for the inter-arm and spiral arm regions 
respectively (Fig. 5). lAllen etltP (|2004 ) find from models 
of photodissociaion regions (PDRs) that CO (1-0) emission 
is detected when the density is 10 cm~^ < n < 10^ cm~^. 
Whilst the spiral arm densities are high enough to be de- 
tected in CO, much of the inter- arm gas in our simulation lies 
below this regi me and w ould be difficult to detect observa- 
tionallv. Kaufman et al.i (1999) also show that the CO to H2 
conversion factor is roughly constant for a 100 cm~^ cloud 
in a low UV field, providing the column density is > 4 x 10^^ 
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Figure 5. The average density is plotted (top) versus azimuth, 
calculated for a ring centred at 7.5 kpc (6 is measured clockwise 
round the disc). The density of the 10^ K gas (dotted), 100 K 
(dashed) (including HI and H2) and H2 (solid) are shown. The 
bottom figure shows the density of molecular hydrogen for indi- 
vidual particles (of similar phases) with time. 
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Figure 6. The typical time over which gas contains over a certain 
fraction of H2 is shown in this figure. For each percentage of H2, 
the average time is calculated from all particles which exceed that 
fraction of H2 during each simulation. The average dynamical 
time for gas exceeding 90% H2 is 5 Myr, increasing to 20 Myr for 
gas exceeding 10% H2. 

cm~^ (comparable with the column density threshold used 
for our clump finding algorithms). Below these densities, or 
with a higher UV field, H2 is significantly underestimated 
from CO observations since CO is dissociated more readily 
than H2. 

Fig. 4 further indicates that the column density of the 
inter- arm molecular gas is typically between 5x 10~^ — 10~^ 
g cm~^. This is an order of magnitude less than the col- 
umn densities of H2 observed for HISA clouds which are 
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Figure 7. The distribution of mass as a function of gas density 
is shown for the cold gas in the inter-arm regions. The points are 
determined after 100 Myr for the single (crosses) and multi-phase 
(squares) simulations. Gas with densities abo ve 10 cm~^ (10~^^ 
g cm~3) is expected to be cold (100-200 K) ( Glover Mac Lowl 
,2006a,). 

asso ciated with CO emiss ion ([Klaassen et alJ l2005l ) . How- 
ever |Klaasse^^t"al? ("2005") also examine HISA clouds where 
there is no CO emission, but infer, by means of radiative 
transfer techniques, column densities of H2 for these clouds 
of around 10~^ g cm~^. We postulate that the hot compo- 
nent of the ISM confines clouds of cold HI to sufficient den- 
sities to maintain a low level of H2. The densities of these 
clouds imply that they will be largely undetected at present, 
but are high enough that the clouds are likely to be cold. 
We note however that the UV flux in these calculations is 
kept constant and we do not include feedback. 

3.4 Resolution effects 

We performed further simulations with 1 million (i.e. 
500,000 cold gas particles for the multi-phase case) and 
250,000 particles to investigate the effects of numerical res- 
olution. The lower resolution results in lower peak gas den- 
sities since with fewer particles the spiral shock is less well 
resolved. Compared with Fig. 3., the peak densities for the 
cold gas decrease by a factor of 3 when 1 million particles 
are used, and a factor of 8 for 250,000 particles for the multi- 
phase simulations. In the single phase simulations, the peak 
densities are similar for 1 million particles compared to 4 
million, but decrease by a factor of 3 for 250,000 particles. 
Consequently the amount of molecular gas decreases with 
resolution, and only a few per cent of the gas is molecular 
when there are 250,000 particles. 

For the multi-phase simulations, there is also a large 
decrease in the inter-arm densities of the gas, the peak den- 
sities shown on Fig. 7 decreasing by a factor of 10 in the 
lower resolution simulations. This reflects that these den- 
sities are dependent on the densities reached in the spiral 
arms and the confinement of this gas by the hot phase as 
it traverses the inter-arm regions. The average H2 densities 
similarly decrease by a factor of 2-3 in the spiral arms com- 
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pared with Fig. 5, and a factor of up to 10 in the inter-arm 
regions. There is httle change in the the inter-arm densities 
for the single phase results (Fig. 7) where most gas leaving 
the spiral arms quickly diffuses to low densities. Apart from 
the change in the peak densities, the overall distribution of 
densities, as shown in Figures 3 and 7, is similar at different 
resolutions. However the peak densities, as well as molecu- 
lar gas densities, represent lower limits both in the arm and 
inter- arm regions. 



3.5 Structure in the cold gas 

We now consider quantitatively the structure in the cold 
gas for the multi-phase medium, and compare this with the 
single-phase cold gas. We compare the properties of clumps 
from each simulation, using 2 clump finding algorithms. 
Both algorithms are grid based, assuming a 2D grid over 
the galactic disc. The first (hereafter CFl) selects clumps 
solely using a density threshold, whilst the second (here- 
after CF2) uses neighbour lists to assign the most dense 
grid cells and their neighbours to clumps. For CFl, the den- 
sity in each grid cell is calculated and all the cells above a 
certain density threshold are selected. Cells that are adja- 
cent are grouped together and particles within these cells 
classed as a clump. CF2 was provid ed by Paul Clark and i s 
based on the clump-find method of IWilliams et al 
The SPH densities are first smoothed onto a 2D grid. The 
grid cells are then ordered by decreasing density. For each 
cell in turn, the cell and its neighbours (subsequently re- 
moved from the list) are identified as a clump, or assigned 
to a previously defined clump. Again a density threshold 
was set which removed low density cells. The properties of 
each clump are then calculated from the grid cells rather 
than the actual SPH particles. Since CF2 selects the most 
dense particles first, density peaks tend to be separated into 
separate clumps where they would otherwise constitute a 
single clump for CFl. In both cases we used the density of 
the molecular gas, which as seen from Fig. 4b largely re- 
flects the underlying structure, and so we can compare the 
properties of these clumps with observed molecular clouds. 

The clumps found using CFl are displayed for the 
multi-phase (Fig. 8) and the single-phase medium (Fig. 9) 
for a subsection of the disc (the same section is shown in 
Fig. 2). Both these figures are taken after 100 Myr. The res- 
olution of the grid for the algorithm is 5 pc, and the density 
threshold for each cell is 100 M©, i.e. equivalent to a surface 
density of 9 x 10~^ g cm~^. Each clump consists of at least 
30 SPH particles (corresponding to a total hydrogen mass of 
5000 M0) which lie within cells above the density threshold. 
Clumps with fewer particles are discarded. We also show the 
corresponding column density plots indicating the molecular 
gas column density (Figs 8 and 9, lower). 

There is a clear difference in structure of this section 
of the disc for the single and multi-phase gas. The multi- 
phase case shows many smaller clumps along the spiral arm, 
whilst the single-phase plot consists of 3 large clumps and 
several smaller clumps. Over the whole disc, there are 8383 
clumps in the multi-phase calculation and 2066 in the single 
phase calculation. The clump finding algorithm also picks 
out several clumps in the inter-arm regions in the multi- 
phase simulation, whereas there are none for the single-phase 
results. These figures further show the general distribution 
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Figure 8. The top panel shows the distribution of clumps found 
using the clump- finding algorithm CFl for a 1.2 by 1.2 kpc sub- 
section of the disc after 100 Myr. The particles corresponding to 
the 100 K gas are plotted in dark blue, the 10^ gas is not shown. 
The particles in each clump are over-plotted in a different colour. 
The lower panel shows the corresponding column density image, 
using the same scale as Fig. 5 with molecular hydrogen overplot- 
ted in red. 



of cold gas. In the multi-phase simulation, the cold gas oc- 
cupies a much smaller filling factor, located in small dense 
clumps, particularly in the inter- arm regions. By contrast, in 
the single-phase case, the cold gas is much more uniform in 
the inter- arm regions. The second algorithm, CF2, generally 
produced smaller clumps, but with similar results. 



3.6 Properties of molecular gas clumps 

We now examine the properties of the clumps from the sin- 
gle and multi-phase simulations. The results shown use CFl, 
but we discuss any differences with the second clump-finding 
algorithm. Again only clumps with > 30 particles are re- 
tained, and the clumps are located after 100 Myr. The mass 
and radii for the clumps from the 2 simulations are shown 
in Fig. 10. At later times the clump masses increase for the 



© 0000 RAS, MNRAS 000.1110 



8 C. L. Dobbs, 1. A. Bonnell 




).6 




Figure 9. The top panel shows the distribution of clumps found 
using the clump-finding algorithm CFl for a 1.2 by 1.2 kpc subsec- 
tion of the disc after 100 Myr. This is taken from the single-phase 
(50 K) simulation. The dark blue points show all the particles in 
the simulation, but the particles in clumps are overplotted in dif- 
ferent colours. The lower panel shows the corresponding column 
density image, using the same scale as Fig. 4. 



imum extent of the x and y coordinates to calculate the 
area. 

Fig. 10 shows the mass versus radii from the clumps 
found in the multi and single phase simulations from a quar- 
ter of the disc. The average smoothing length over all the 
particles is 30 pc, but reduces to 10 pc for particles above the 
cutoff density used for the clump-finding algorithm. There 
are far fewer clumps in the single-phase simulation - less 
than 1/3 of the multi-phase simulation. However there are 
comparatively more larger clouds, as also indicated by Figs 8 
and 9. The total mass of the clumps in Fig. 10 for the single- 
phase simulation is 5.7 x 10^ compared to 1.4 x 10^ 
M0 for the multi-phase simulation, over a quarter of the 
disc. The most massive clump after 100 Myr in both simu- 
lations was ^ 3 X 10^ M©, with a radius of 150 pc. For the 
multi-phase simulation, the fraction of molecular hydrogen 
increases with time, so after 200 Myr, the maximum cloud 
mass is '-^ 6 X 10^ M©. 

The clumps found using the second algorithm, CF2, fol- 
low a similar distribution to that shown on Fig. 10. How- 
ever there is less scatter and fewer very massive clouds. The 
clumps from both simulations, and from both algorithms 
follow a M oc dependence, implying the clouds have con- 
stant surface densities of ~ 8 M© pc~^ (2 x 10~^ g cm~^). 
For comparison, the density criterion of 4 M© pc~^ is also 
shown. The apparent constant surface density of the clumps 
is a consequence of the density threshold applied to find the 
clumps. Increasing the threshold 10 times moves the dis- 
tribution of clumps to correspondingly smaller radii. The 
surface de nsities of molecular clouds from observatio nal re- 
sults (e.g. iHever k TerebevI Il998l : iBlitz et all l2007l ) could 
similarly represent a selection effect based on the threshold 
antenna radiation temperature, T^. 

Since we do not include any feedback effects, most of 
the clumps are situated in the mid-plane of the disk, with 
average scale heights of < 10 pc. However approximately 
10% of the clouds, generally those of lower mass, have a 
centre of mass between 50 and 100 pc above or below the 
plane of the disk. 

The resolution of most of the clumps was insufficient to 
calculate a velocity dispersion. However the better resolved 
clumps (with > 100 particles) contain velocity dispersions 
of a few km s~^. We also calculated the virial parameter. 
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multi-phase simulation increase, as there is a higher fraction 
of H2, while those of the single phase simulation are similar. 

For both clump-finding algorithms, we take a column 
density threshold 9 x 10~^ g cm~^. For CFl, this corre- 
sponded to taking a grid of resolution 5 pc and requiring 
that each cell contains at least 100 M© (of H2). For CF2, 
we apply the same grid resolution, and a minimum column 
density of 0.001 g cm~^ of H2. With CFl, each clump cor- 
responds to a group of SPH particles. The radii shown in 
Fig. 10 are determined by calculating the radius which con- 
tains 3/4 of the total clump mass, and thus the mass plotted 
is 3/4 of the total clump mass assigned by the particles. Us- 
ing the full or 1/2 clump radii and masses merely shifts the 
distribution of clumps to higher or lower masses and radii. 
For CF2, we use the effective radius from the grid cells which 
make up each clump, i.e. Veff = Areaji^ using the max- 



for the clumps. In both simulations, most of the clumps are 
unbound. However the more massive clumps tend to be more 
bound, with cxq typically < 10. 

Fig. 11 shows the mass spectrum from the 2 simulations, 
for the total number of clouds found over the whole of the 
disc (of > 30 particles) with CFl. Assuming a mass spec- 
trum dN/dM oc ^ we find for the multi-phase simula- 
tion, a — 2. 35 ±0.15. The clumps from the single-phase sim- 
ulation produce a much shallower spectrum, ofa = 1.6=b0.3. 
At later times, the spectrum for the single-phase simulation 
is similar, although the spectrum for the multi-phase simu- 
lation becomes shallower, with a ^ 2.15 after 200 Myr, since 
there are more larger clouds at later times. Using CF2, the 
larger clumps divide into smaller clumps which truncates 
the clump masses at 1.5 x 10^ M©. Consequently the value 
of a increases to a ^ 3.3 and a ^ 2.6 for the multi and 
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Figure 10. The mass and radius of the clumps are plotted us- 
ing CFl over a quarter of the disc. Lower mass clumps continue 
this distribution, but are not shown. The clumps are selected at 
100 Myr. The dashed hue indicates E = 9 x 10"^ g cm^^ (4 Mq 
pc-2). 



single-phase simulations respectively. The form of the mass 
spectrum thereby appears to depend on the details of the 
clump-finding algorithm. What can be concluded however 
is that the multi-phase medium produces a steeper mass 
spectrum and fewer high mass clouds. Observations indicate 
that a lies in the range a = 1.5 — 1.8 (jSolomon et al.lll987l : 
iHever Terebevlll998[L for the Galax y and up to a = 2.5 
for external galaxies ([Blitz et"al]l2007l ). 

Fig. 12 shows the clump mass spectra for the multi- 
phase simulation with 6 million, 1 million and 250,000 parti- 
cles. The clumps are again found using CFl, using the same 
column density threshold (hence the clumps lie within the 
distribution in Fig. 10), but increasing the grid-size for the 
clump-finding algorithm. All the clumps contain > 30 parti- 
cles. The slope becomes shallower as the resolution decreases 
and with 1 million particles the total amount of H2 is sig- 
nificantly reduced to 10%. At lower resolutions, less massive 
mass clouds cannot be detected, although the upper mass is 
not significantly reduced compared to the higher resolution 
simulations. For single phase calculations the slope of the 
spectrum was the same for 4 and 1 million particles, but 
shallower for 250,000 particles. In summary, we stress that 
the results presented here should be taken as lower limits 
on the fraction of molecular gas present, both in the spiral 
arms and in the inter arm regions. 



3.6.1 Inter-arm and spiral arm molecular clouds 

As mentioned previously, the multi-phase simulation pro- 
duces a much greater degree of inter- arm molecular gas and 
a few of the clumps in Fig. 8 are located in the inter-arm re- 
gions. We show the distribution of molecular clouds over the 
whole of the disc, using CFl, in Fig. 13. There are a signifi- 
cant number of clouds in the inter-arm regions, especially in 
the inner regions of the disc where the spiral arms are closer 
together. For the outer parts of the disc, the inter-arm clouds 
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Figure 11. The clump mass spectrum is shown for clumps from 
the multi-phase simulation (solid) and the 50 K simulation (dot- 
ted). Both are determined after 100 Myr and use CFl. The masses 
indicated are the mass of H2. The imposed limit of 30 parti- 
cles corresponds to a total (HI and H2) mass of approximately 
10^-^ M0. Allowing that the actual mass of H2 is a fraction 
of the total, the resolution limit corresponds with the peak of 
the distributions, below which the sample is incomplete and the 
number of clumps rapidly drops off. The slope corresponds to 
dN/dM oc M~^-^^ for clumps from the multi-phase simulation 
compared to dN/dM oc M~^'^ for the single phase results. 



tend to be situated on the edge of spiral arms. We estimate 
the ratio of inter- arm to arm clouds is approximately 1:7 
(decreasing to 1:10 for the 1 million particle simulation, and 
no inter-arm clouds at lower resolution). The distribution of 
clouds and the inter-arm ratio is very similar for the second 
clump-finding algorithm. By contrast, for the single phase 
simulation, the clouds all lie along the spiral arms. Con- 
sequently observational- sty le plots of the molecular clouds 
over a range of longitudes fPobbs et al.ll2006l Fig. 14 k 15) 
show a greater degree of scatter in the inter-arm regions. 

We also considered the number of more massive clouds 
in the inter-arm and arm regions. Approximately 3% of 
inter-arm clouds were > 10^ Mq compared to 6% of clouds 
in the spiral arms. The cloud mass spectra for clouds lo- 
cated in the spiral arms and inter- arm regions are displayed 
in Fig. 14. The spectrum for the spiral arm clouds is very 
similar to the spectrum for the total number of clouds, al- 
though slightly shallower with a = 2.25 =b 0.15. The spec- 
trum for the inter-arm clouds however is much steeper with 
a = 2.9 =b 0.4, and a cut off in mass at around 3 x 10^ M©. 

The number of inter-arm clouds is sensitive to the pho- 
todissociation rate, so these figures only provide a rough 
indication. The photodissociation rate may be lower in the 
inter-arm regions, increasing the number of inter-arm clouds. 
On the other hand, the photodissociation rate will be higher 
where massive star formation occurs, which with the effects 
of feedback may disrupt molecular gas clouds before they 
enter the inter- arm regions. 
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Figure 12. The clump mass spectrum is shown for the multi- 
phase simulation using different resolutions. The clump mass is 
the mass of H2. The total number of particles is indicated for each 
simulation. 



Figure 14. The clump mass spectrum is shown for the clumps 
in the spiral arms (dotted) and inter-arm regions (solid) from 
the multi-phase simulation. The clumps are selected at 100 Myr, 
and again the minimum (total) mass resolution of the clumps is 
103-5 Mo. 




Figure 13. The distribution of clumps is displayed for the multi- 
phase simulation, using CFl after 100 Myr. 



4 CONCLUSIONS 

We have performed hydrodynamic simulations of a multi- 
phase medium subject to a spiral potential. We find that in- 
cluding hot gas in addition to cold leads to higher densities 
in the cold gas. Consequently the amount of molecular gas 
increases in the disc, with 1/4 of the cold gas molecular after 
100 Myr, approximately double that of corresponding sin- 
gle phase calculations. This fraction further increases with 
time. The increase in H2 is most striking in the inter-arm 
regions, with molecular clouds surviving into the inter-arm 
regions without being destroyed by photodissociation, al- 
though both feedback, and heating and cooling of the gas 
are neglected. For the whole disc, approximately 1/7 of the 



molecular clouds lie in the inter-arm regions, but this ratio 
decreases with radius as the separation of the spiral arms 
increases. The spiral arm H2 densities increase to values 
comparable with average GMC densities, whilst the typi- 
cal inter- arm densities are below those detected with CO. 
We therefore propose that a multi-phase medium predicts 
a population of cold clouds of HI and molecular hydrogen, 
which are currently undetected. 

The multi-phase simulations also show more structure 
in the cold gas, with smaller clumps present compared to the 
single phase calculations. Consequently the mass spectrum is 
steeper for the multi-phase clouds, although the value of the 
exponent is found to depend on the clump-finding method. 
The maximum cloud mass is 6 x 10^ i n the multi-phase 
simu lation, comparable to GMCs in M33 (|Engargiola et aP 
l2003l \ but an order of magnitude lower than the most mas- 
sive GMCs observed in the Milky Way. The free-fall times 
for the molecular clouds are esimated to range from 5 Myr, 
when over 90% of the gas is molecular, to 20 Myr when over 
10% of the gas is molecular. The lifetimes of the clouds are 
estimated to be 20 Myr for when over 90% of the gas is 
molecular, and 80 Myr when over 10% of the gas is molecu- 
lar, both of which correspond to roughly 4 free-fall times (5 
Myr for > 90% and 20 Myr for > 10% H2). 

The main factor which will affect these results is the 
inclusion of feedback, particularly in determining how much 
molecular gas is retained in the ISM following star forma- 
tion. This may lead to less inter-arm molecular gas, but 
would also introduce the possibility of triggered molecular 
cloud formation in the inter-arm regions. Furthermore, us- 
ing lower scale heights to determine the photodissociation 
rate could also reduce the amount of H2.. Thirdly we do not 
include heating and cooling of the different phases. Rather 
we have assumed a fixed composition of the ISM to approxi- 
mately agree with observations. Previous simulations of the 
ISM with heating and cooling sh ow that a cold phase of gas 
is sustained, wit h 20-30% of gas (lAudit Hennebellell2005l) 
to 60% or more (jWada Normanll200ll : IPiontek Ostriked 
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I2OO5I ) in the cold regime. We leave a more complete treat- 
ment of the thermodynamics of the ISM and spiral shocks 
for future work. 
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